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Abstract. Seismology of single S Scuti stars has mainly been inhibited by failing to detect 
many of the theoretically predicted pulsation modes, resulting in difficulties with mode 
identification. Theoretical and observational advances have, however, helped to overcome this 
problem, but the following questions then remain: do we know enough about the star to either 
use the (few) identified mode(s) to probe the structure of the star? or improve the determination 
of the stellar parameters? It is now generally accepted that for the observed frequencies to be 
used successfully as seismic probes for these objects, we need to concentrate on stars where we 
can constrain the number of free parameters in the problem, such as in binary systems or open 
clusters. The work presented here, investigates how much is gained in our understanding of the 
star, by comparing the information we obtain from a single star with that of an eclipsing binary 
system. Singular Value Decomposition is the technique used to explore the precision we expect 
in terms of stellar parameters (such as mass, age and chemical composition) as well as how these 
parameter uncertainties propagate to the Luminosity- Temperature (L-T) diagram. This work 
shows that the information content of the binary system provides sufficient constraints on the 
models so that the mode can be used to probe the star's structure. 



1. Introduction 

d Scuti stars are a class of pulsating stars located on the Hertzsprung-Russell Diagram on or 
around the Main Sequence and intersecting the instability strip. They are 1.5 - 2.5 Mq stars 
pulsating often in one main dominant oscillation mode or many lower amplitude pulsation modes. 
They have been interesting targets seismologically, because the oscillation amplitudes often reach 
tenths of magnitudes, and we understand their stellar structure relatively well, so seismology can 
allow us to probe details of the microphysics such as energy transport mechanisms, convective 
core overshoot, as well as other less well-developed theories such as rapid rotation. 

The main setbacks that 5 Scuti seismology face (as well as other pulsating stars) are 1) the 
fundamental stellar parameters are not well-enough constrained to allow the few pulsation modes 
to probe the structure, and 2) rapid rotation causes each of the degrees I to split into 21 + 1 m- 
modes, making mode-identification a difficult task pp. Indeed these problems are not exclusive 
nor exhaustive. This has prompted authors to look towards objects where at least one of these 
problems can be eliminated [21 [31 HI [5]. Observing stars where the number of free parameters 
is constrained, such as in open clusters or multiple systems is a possibility for overcoming these 
obstacles. 



In order to use seismology to probe the interior of a star, the parameters of the star need 
to be known quite well, for example, the mass should be known to 1-2% [6]. The observables 
from a binary system provide strict constraints on the parameters of the component stars. If 
the binary is an eclipsing and spectroscopic system, the absolute values of the masses and radii 
can be extracted to 1-2% (e.g. [3 IB])- 

The objective of this study is to find out if the uncertainties in the stellar parameters can 
be reduced, so that seismology can be applied to those stars that exhibit one or few pulsation 
modes. We look at the particular case of a pulsating star in an eclipsing binary system and 
compare the parameter uncertainties with those of an isolated star. This study quantifies how 
well the stellar parameters can be extracted in various hypothetical systems. 



2. Methods 

The mathematical basis of this study lies in the application of Singular Value Decomposition 
(SVD) techniques to physical models. This technique has been used in previous studies, such as 
Ell E] and it is being applied to many other areas of astrophysics and science, because of its 
powerful diagnostic properties. Intricate details of this mathematical technique are elaborated 
upon in the above mentioned publications as well as in [11] . Here we give the basic equations 
to enable an understanding of this work. 



2.1. Singular Value Decomposition 

SVD is the decomposition of any M x N matrix D into 3 components U, V T and W given by 
D = UWV T . V T is the transpose of V which is an N x N orthogonal matrix that contains the 
input basis vectors for D, or the vectors associated with the parameter space. U is an M x N 
orthogonal matrix that contains the output basis vectors for D, or the vectors associated with 
the observable space. W is a diagonal matrix that contains the singular values of D. 

The key element to our work is the description of the matrix D. Here we define D to be a 
matrix whose elements consist of the partial derivatives of each of the observables with respect 
to each of the parameters of the system, in function of the expected measurement errors on each 
of the observables: 

Here B{ are each of the i = 1, 2, ...M observables of the system, with measurement or expected 
errors 6j, and Pj are each of the j = 1, 2, N free parameters of the system (see section [2~2l for 
discussion on the observables and the parameters). 

By writing the design matrix in function of the measurement errors, we provide a quantitative 
description of the information content of each of the observables for determining the stellar 
parameters and their uncertainties. 

Supposing that we are looking for the true solution Pr of the system. By starting from 
an initial close guess of the solution Pq, SVD can be used as an inversion technique by 
calculating a set of parameter corrections <5P that minimizes some goodness-of-fit function: 
5P = VW _1 U T £B, where <5B are the differences between the set of actual observations O and 
the calculated observables Bo given the initial parameters Po • W is a modification of the matrix 
W such that the inverse of the values below a certain threshold are set to 0. The formal errors 
are comprised of the sum of all of the ~Vk/ w ki where each Vk/ W k describes the direction and 
magnitude to move each parameter, so that the true solution Pr and formal uncertainties can 
be given by 

Pr = Po + vw- 1 ^ f ±Xi ± ^ ± ... ± x^y (2) 



The covariance matrix C consequently comes in a very neat and compact form: 

= £ Wfc (3) 

and the square roots of the diagonal elements of the covariance matrix are the theoretical 
parameter uncertainties: 




2.2. Observables, Parameters & Models 

We describe a single 5 Scuti star by a set of parameters or ingredients. The ingredients for the 
stellar model are mass M, age r, rotational velocity v, initial Hydrogen and heavy metals mass 
fraction X and Z where X + Y + Z = 1 and Y is initial Helium mass fraction, and mixing-length 
parameter a where applicabl^]. The distance to the object d is also included as a parameter. 

For a binary system, the additional parameters are the system properties: separation of 
components a, eccentricity of orbit e, longitude of periastron u and inclination of orbit i. 
Fortunately, both stellar components in a binary system share the parameters r, X and Z, 
so then the individual stars differ mainly by M and v. The parameters of the binary system of 
this study are given in Table [TJ 

The observables are the measurable quantities of the system. These include things such as 
radius, effective temperature, gravity, metallicity and parallax for a single star system. For a 
binary system, the observables include effective temperature ratio, relative radii, radial velocities 
and orbital period. The Aarhus STellar Evolution Code [12j is used to calculate the stellar 
evolution models. This code uses the stellar parameters as the input ingredients, and returns a 
set of global stellar properties such as radius and effective temperature, as well as the interior 
profiles of the star such as mass, density and pressure. Oscillation frequencies for this rapidly 
rotating star are calculated using MagRot [131 E] • Using the global stellar properties and the 
distance to the star, SDSS [15] magnitudes and colours in various filters are evaluated using the 
Basel model atmospheres |16| . Most of the binary observables are calculated analytically using 
the combined stellar and binary parameters. In this work we refer to the non-seismic data as 
the classical observables. Most of the classical observables of the binary system are given in 
Table E 

A clear distinction should be made between the input parameters of the system and output 
measurable quantities, the observables. So to discriminate between the errors in both parameters 
and observables, we shall denote the (derived) parameter uncertainties by a, while e is reserved 
for the observable errors. 

Both luminosity and effective temperature can be observables, but in Section [3. 21 the derived 
uncertainties of these properties are discussed. This refers specifically to the calculated error 
boxes in the luminosity-temperature (L-T) diagram, and here their uncertainties will also be 
denoted by a. 

3. Results 

The theoretical uncertainties in each of the parameters of Ma, t, X and Z are calculated as 
a function of error in radius and as a function of error in colour e^_ z ), using Equation 2] 
coupled with the observable errors given in Table [TJ Consequently the theoretical uncertainties 
in luminosity and effective temperature are calculated using Equation [21 for three different values 



For masses larger than about 1.5 Mq the outer convective layer is relatively thin, so the observables of the star 
are not very sensitive to the value of the mixing-length parameter. 
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of and three different values of effective temperature error £T cff . The results are shown in 
Figures [T] and [5] below. 

3.1. Parameter Uncertainties 

Figure [1] shows the theoretical uncertainties (a) in Ma, t, X and Z of a pulsating star as 
a function of error in the radius (left panels) and as a function of error in the photometric 
colours (right panels). Note that the left panels do not include photometric information, i.e. 
no colours nor magnitudes. The dashed lines show the results for a single star (the observables 
are radius, effective temperature, gravity, metallicity, ...) while the solid lines show the results 
for a component of a binary system (observables are those of the single star and the binary 
observables). The lines with the diamonds include one identified mode as well as the classical 
observables. The results for the other stellar parameters are not shown, because the four 
aforementioned parameters and va are responsible for determining the model structure of the 
pulsating star. <j(va) is is usually independent of and e^_ z y, it is determined mainly by the 
observables vsini and i from a combination of spectroscopy and the photometric light curve. 

For the single star without an identified mode (dashed lines no diamonds), the parameter 
uncertainties remain at a large constant value as a function of radius (left panels) but do decrease 
slightly with improved photometric data (right panels). Only when seismic data are included 
for the single star system (dashed lines, diamonds) and the observable errors are small, the 
parameters are constrained to a usable amount. 

By comparing the solid lines with and without diamonds in Figure HJ it can be seen that 
the addition of one identified mode makes almost no difference to the parameter uncertainties 
for the binary system. This implies that there is enough information provided by the binary 
constraints to sufficiently determine the stellar parameters. In this sense, the identified mode is 
redundant information, and thus can be used maybe to test the interior of the star. 

Including photometric information (right panels) provides an interesting result: the 
information provided by the single star system can supersede that of the binary system for 
r and Z . This is because the colours are uncontaminated by a component star. This only 
happens at very small measurement errors, and only when an identified mode is included for 
the single star. 




Figure 1. Theoretical uncertainties (a) in mass M, age r, initial hydrogen X and metal Z 
content as a function of observable error. The left panel shows the uncertainties as a function 
of observable radius error, here no photometric information has been included. The right panel 
shows the uncertainties as a function of observable error in photometric colour. The dashed and 
continuous lines show the results for the single star and the binary system respectively, those 
with diamonds show the results when an identified mode is included in the set of observables. 




Figure 2. The theoretical error boxes for luminosity and effective temperature. The dashed 
lines represent the results for the single star while the continous lines represent the results for 
the binary system. The left panel shows the results while reducing the error in radius, and the 
right panel shows the results while reducing the error in effective temperature. 



3.2. Luminosity- Temperature Error Box 

The correlation matrices come in a compact form when using SVD. This then allows a calculation 
of the theoretical uncertainties in both effective temperature and luminosity (L-T) (Equation 
Figure [2] shows the theoretical error box in effective temperature and luminosity for a single 
star system [dashed lines) and a binary system {solid lines). The observables do not include 
photometric information, and for the single star an identified mode is included^!, while for the 
binary system no seismic data is included. 

3.2.1. Single Star Observe how the error box reduces significantly while reducing the error 
in the radius observable (left panel). The uncertainty in T e g also reduces slightly. The right 
panel also shows that by reducing the errors in the observable T e g, an expected corresponding 
reduction in the uncertainties in T e g is noted. The €t cH of 200, 100, and 50 K, produces a cr(T c g) 
of 250, 110, and 50 K. The fact that these uncertainties are reproduced also gives confidence in 
this method. cr(L+) changes slightly as a function of er efi , its value is determined mostly by the 
error in the radius observable (2%). Looking back to the left panel, we see that interpolating 
between 1% and 3% ej^ produces a cr(L+) = 0.5 L© for = 2%. This is the value that is 
shown in the right panel. 

3.2.2. Binary System For the binary system (solid lines), no identified mode is included. The 
error box for the binary system does not reduce while reducing the errors in the radius, because 
of the small uncertainties in these parameters. However, the error box does reduce when the 
error in effective temperature is reduced, reproducing accurately the input 6t cS of a(T e ^) = 200, 
100, and 50 K. o~(L+) does not decrease in either panel, because the mass is well-determined for 
the binary system and provides this narrow constraint on L+. 

In all cases, note that the constraints provided by the binary system without an identified 
mode are more effective than those from the single star when an identified mode is included. 

2 Figure [1] shows that the parameters are not constrained for the single star if the identified mode is not included. 



4. Conclusions 

This study investigated whether the uncertainties in the stellar parameters of a pulsating 
component in an eclipsing binary system were sufficient so that an observed pulsation mode 
could be used to test the physics of the stellar interiors. Additionally we studied the information 
content of a pulsating star in a single star system to quantify how much is gained in terms of 
precision in parameters and size of the L-T error box by observing the star in a detached eclipsing 
binary system. The conclusions are summarized as follows: 

• A single star system without an identified mode remains poorly understood when 
observables such as the radius or colours are poorly measured. The parameter uncertainties 
are too large to correctly place the star in the L-T diagram. 

• A binary system without seismic information provides better constraints than the single 
star system when an oscillation mode has been identified. 

• Reducing the size of some observable errors has little or no impact on the parameter 
determinations for the binary system, because these parameters are already well 
constrained. 

• The tight constraints provided by the binary system for the stellar parameters reduces the 
size of the error box in the L-T diagram significantly. 

• By carefully constraining the parameters of the star, just as an eclipsing binary system 
allows us to do, an accurate estimate of the stellar model under study can be obtained. 
This allows the redundant observables (like an oscillation mode) to be used exclusively to 
test the physics of the interior of a star. 
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